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ON THE TIDAL EFFECTS IN THE MOTION 
OF ARTIFICIAL SATELLITES 

Peter Musen 
Ronald Estes 

Goddard Space Flight Center, Greenbelt, Maryland 
ABSTRACT 

The general perturbations in the elliptic and vectorial elements of a satel- 
lite as caused by the tidal deformations of the non-spherical Earth are developed 
into trigonometric series in the standard ecliptical arguments of Hill-Brown 
lunar theory and in the equatorial elements oj and Q of the satellite. 

The integration of the differential equations for variation of elements of the 
satellite in this theory is easy because all arguments are linear or nearly linear 
in time. 

The trigonometrical expansion permits a judgment about the relative sig- 
nificance of the amplitudes and periods of different tidal "waves" over a long 
period of time. 

Graphs are presented of the tidal perturbations in the elliptic elements of 
the BE-C satellite which illustrate long term periodic behavior. The tidal effects 
are clearly noticeable in the observations and their comparison with the theory 
permits improvement of the "global" Love numbers for the Earth. 
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BASIC NOTATIONS 

V ,F, D and F - the arguments of the lunar theory, 

G - the gravitational constant, 

M - the mass of the Earth, 
m' - the mass of the Moon, 
m" - the mass of the Sun, 

7 - the geocentric position vector of the satellite. 


7° - the unit vector in the direction of r , 

7' - the geocentric position vector of the Moon, 
r' = |7'1 , . 

7'“ - the unit vector in the direction of 7' , 

7" - the geocentric position vector of the Sun, 

_« - |-;r» 1 

r Ml, 

7"® - the unit vector in the direction of 7" , 

a' - the mean distance of the Moon from the Earth. 
Its numerical value is defined in such a manner 
that the constant part in the expansion of the 
parallax a' /r' is equal to one. 

a" - the mean geocentric distance of the Sun, 

S' - the angle between 7 and 7' , 

S" - the angle between 7 and 7" , 

K, /X, u, k' , /x', i/', k", jj." v" - the components of 7^, 7'“ and 7*® respectively 

in the equatorial system, 

e - the eccentricity of the orbit of the satellite, 

i - the orbital inclination of the satellite. 
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In the beginning of the exposition these symbols designate the osculating 
elements. Toward the end they designate the mean elements. 

n - the longitude of the ascending node of the satellite, 

o) - the argument of perigee of the satellite, 

77 - the longitude of the perigee of the satellite, cj + Q, 
c = cos i/2, 
s = sin i/2, 

SM - the perturbations in the mean anomaly of the 
satellite caused by tides, 

SO - the perturbations in 0 caused by tides, 

S 77 - the perturbations in 77,. 

Si- the perturbations in S i caused by tides, 

n - the mean motion of the satellite, 

a - the semi-major axis of the satellite's orbit, 

R - unit vector normal to the orbital plane of the 
satellite, 

Uj, U2, U3 - components of R in the equatorial system, 

Suj, Suj, SU3 - perturbations in Uj ,Uj ,03 respectively caused 
by tides, 

6 - the obliquity of the ecliptic, 
e - the eccentricity of the Earth’s meridian, 
kj, k^, • ' ‘ - Love numbers, 

R - the equatorial radius of the Earth, 
a = R/a 
a' = R/a’ 
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INTRODUCTION 


The precision of observations of artificial satellites has now reached the 
stage where consideration must be given to the effects in their motion as caused 
by the tidal deformations of the Earth. The comparison of tidal theory with ob- 
servations will permit us to determine the '’global” Love numbers for the Earth. 
They appear as coefficients in the expansion of the tidal potential and describe 
the average elastic properties of the Earth as "seen” from outer space. In 
connection with this topic the works by (Kozai, 1965), (Fisher and Felsentreger, 
1966) and (Newton, 1967) must be mentioned. The tidal effects in the motion of 
satellites are very small, a few seconds of arc, and their computation, as well 
as the determination of Love numbers, requires a great care in the expansion 
of the tidal disturbing potential and In the numerical computations. 

We discard the usual trigonometric expansion of the tidal potential with the 
angular equatorial elements of the Moon as arguments. This typo of expansion 
is very simple in form, but the analytical integration of the resulting differential 
equations is difficult because the mean equatorial node of the lunar orbit does 
not change linearly with time, but oscillates between two fixed limits. 

We suggest instead the expansion of perturbations into trigonometric series 
with the ecliptical arguments 't' , F, D andP of the lunar theory and with the 
equatorial elements oj andfl of the satellite. 

The arguments -t, , F, D and P are very nearly linear with respect to 

time and, as a consequence, the integration becomes much simpler and will be 
valid for a longer interval of time. 
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We expand the tidal potential into a sum of products of two harmonic functions, 
the first depending on the position of the satellite and the second on the position 
of the Moon (Sun). 

In order to facilitate the programming we represent the harmonic functions 
in Maxwellian form as polynomials in equatorial components of unit vectors. 

We obtain the trigonometric expansions of the functions , ij.\ v' and 
(a*/r')^ which appear in the tidal potential by making use of the standard devel- 
opments of the lunar parallax a'/r', latitude 0^ and related to the lunar 
longitude by 

= + D+ r + 8X^ 

into Fourier series (Brown, 1919), (Eckert, 1966), (Eckert et al., 1969), (Siry, 
1970), (Amer. Eph., 1972). The accuracy of these series presently is sufficient 
to treat the tidal effects even for a relatively long period of time. 

The equatorial components of the lunar unit vector are 

k‘ = cos cos 

/x' = s in cos cos e - s in /3^ s in e 

v' = s in cos e + s in cos s in e 

where e is the obliquity of the ecliptic. Expressions for sin cos , sin 
and cos are obtained from 

sin 5X, = hk^ - i + • • • 

cos = 1 -1 (8X^)2 + ± 
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cos = cos ('t' + D + P) cos - s in ('t* + D + P) sin 8>^ 

s in = sin ('t' + D + P) cos 8A^ + cos (-t' + D + P) sin 8A^ 

sin /3, = + • • • 

cos /3, = 1 - y/3j + • • • . 

It is sufficient to retain in 8X^, /3^ and a'/r' only the terms with amplitudes 
larger than 0.5 x so that the resulting accuracy of the expansion of the 
harmonic functions depending on the Moon's position will be 1 x 10~^. The ac> 
curacy of the expansion of the tidal potential is then 6 x 10~^^ and the formal 
accuracy of perturbations is about 6 x 10"* or OllOOl, 

The expansion of the tidal potential is obtained on the machine. PL/I was 
used by one of us (R. E.) to program the operations with Fourier series. 

The Fourier series depending on the moon's position have purely numerical 
coefficients and the decision concerning the relative in^>ortance of terms is 
handled by programming logic. All terms with amplitudes too small to be observed 
are rejected automatically. The machine calculations, in addition to standard 
checks for numerical accuracy, are subject to other chocks on procedure and 
accuracy by the formal properties of the e}q)ressions. 

Only the long period tidal effects can be easily observed and, consequently, 
we suppress all the terms with periods nearly equal to or less than the period 
of revolution of the satellite by averaging the force function over the satellite's 
orbit. The coefficients of the averaged expansion are simple rational functions 
in e, /l - e^, cos i/2, and sin 1/2. 
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To analytically continue the tidal potential into extraterrestrial space we 
need a suitable approximation to the Earth's surface. 

Rather then assumming sphericity, as is usually done, we select an ellipsoid 
of rotation as this approximation. It is of interest to note that oblateness will 
produce near resonance effects in the neighborhood of some critical angles, 
like 46?4, 58?1, 63f4, 69.°0 and 73.°1. 

In observations these resonances cannot be separated from the similar 
resonances caused by the direct lunar attraction. In addition, they are hardly 
noticeable when the eccentricity is small. It is also of some Interest to note 
that the oblateness causes a term of the form 



to appear in the extraterrestrial tidal potential, where e designates the eccen- 
tricity of the meridian and T is a harmonic function of the Moon's position. 

The term is very small, but it can gain importance as the requirements put 
on the accuracy of the determination of GM gra dually increase. 

The aggregate of terms associated with k 2 and not containing rr = cu -fQ in 
the arguments constitutes the largest and the most important portion of the ex- 
terior tidal potential. It defines the so-called "main problem" in the theory of 
tidal effects. Obviously the main problem in the present exposition is not associated 

with the spherical approximation to the Earth surface, as the oblateness enters into 
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the coefficients of the expansion. The semi-major axis and the eccentricity 
remain invariant under the influence of the tidal forces defined by the main 
problem (Kozai, 1965). 

The lag of tides in time can be determined from the shift in phase of the 
equatorial node of the lunar orbit (Kozai, 1965). 

However, in the process of expansion into Fourier series this shift can be 
associated with the element fi, because n' always appears in combination Q - Q' . 
A similar shift exists in the arguments F and D. 


THE EXTERIOR TIDAL POTENTIAL 

The expansion of the lunar tidal force function on the surface of the Earth is: 


V' = Gi 


i' k, J_ cos^ S' - 1 ^ 

' r'3 \2 2 1 

Gm' k, — — cos’ S' - — cos S'^ 

^ r"» \2 2 I 


( 1 ) 


Similarly for the Sun we have: 


V = Gm^k, 


+ Gm'k. 


/ 3 2 c* ^ \ 

, x COS'* S - - 

r«3 \2 2 / 


— f- 

r-4 \2 


cos’ S* - COS S* ) + 


( 1 ’) 


Making use of the expansion of the disturbing function for the direct pertur- 
bations in the motion of an artificial satellite into a series of products of spherical 
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harmonics (Musen, Bailie and Upton, 1961) we can re-write (1) in the form: 


V' = 


Cm* 


^20 


(®21 ^21 ^21 ) 


21 ‘'21- 


(aj2 Cj2 


^22 ^22 


)) 


Cm' 


,*4 


Uaa, C;, 


^31 SJi) 


^32 Sia 


+ (^33 ^33 ^33 ^33^ * (®34 ^34 + ^34 ^34^^ + 


where we let 


®20 = + |r2(l - 3;/2) 

®21 = + b2j=|r2\/x, 


( 2 ) 


and 


= + 3 


bj2 = + 3 \ i/, 


aji = b3, = +|rV(l -5^2). 


bjj = + i z/ (3 - 5 ^ 


^33 = +1 r 2\(\2 _ 3 ^ 2 )^ b ^3 = + I r 2 / i ( 3 A 2 _ ^ 2 y 


15 3, 

^34 = 


^34 = + ^ r2i.(\2 _ ^2), 




(1 -3u'^) 
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The extensions A. . , Bj^ of the harmonic functions a , b j ^ into extraterrestrial 
space are respectively 


. 1 2 L 55 2\ R® A 2 R^ * 

= T * i‘-42^) 

^200 " 4 (y ^ (1-3 v^) , 

*201 (f y (3-30^2 *35 

*21 = (‘ - I -'') ^^ l 0 + ''-^*211 

*210 = +J (7] 
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^31 " ■*■” 7 ^ 310 ’ 




^310 = +^ (- ) 


A 32 - 0 , 


®32 “ ■^— 7®320’ 


■^320 


- m '’ 


(3- 51-2), 


, . R’ . 

"33 ~ ~~7 ”330 > 

a^ 


^330 = (7') ^(^^- 3 /x 2 ) 


^33= +^«33l 

a* 


B 331 =4 


8 Vr 


. _ R\ 

^34 “ ■'■“7 ”340’ 
a^ 


A3,o=4(7y^M.. 


B34 = +__ B 34 q , 

a^ 


B ,40 = + ^ V (\2 _ ^ 2 ) 


'340 
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The correctness of these extensions can be checked directly. 

Replacing ^ in (2) by the corresponding A ^ . and taking 

= GM 

into consideraticm we obtain the following decomposition of the tidal potential 
acting on the satellite: 

V' 

tn^ e» k, WJ, 

I 

+ n* a ’ kj W30 , 

where 

^20 “ ( ^ ^ ^00 ^0 

+ ^1 - (^210 ^21 ®210 ^ 21 ^ 

+ (^20 ^22 ■'■ ®220 ^ 22 ^ ’ 
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^^21 " ^201 ^20 + ^21 ® 2 ll ^ 21 ^ 

+ (A221 C22 + B 22 I Sjj)* 

^30 ~ (^310 ^31 ■*■ ®310 ^ 31 ^ ■*■ ®320 ^32 

+ (^330 ^33 ■•■ ®330 ^ 33 ^ ■*■ ^^340 ^34 "*■ ®340 ^ 34 ^ ' 

It is of interest to note the presence of the factor GM/r in the first term. 
The expansion of the factor associated with GM/r consists of a constant term 
and a set of long period terms depending upon the arguments -t' , F, D and P 
of the lunar theory. Consequently, the contribution of the oblateness to the tidal 
potential is a term which may influence tiie accurate determination of the basic 
quantity GM. 


The main part of V associated with the "main problem" of the determination 
of tidal effects is, of course, Because of the additional parallactic factor 

in W'q we have used only the spherical approximation to the Earth's surface to 
obtain the extension of terms depending upon into extra-terrestrial space. 

We have: 


1 


[ 4 \ 

42 j 


+ 1 C' 

^ 2 V 14 j 210 

+ ^ "if sc(c2 -s2)C'2o 


(1 -e2)'3/2, 
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where 


C' = C' 

''200 '' 20 ’ 

Cjjo = Cji COS 2 n + Sjj sin 2 Q, 

C '20 = C'j cos n - S 22 sin 0, 

are the trigonometrical series in^t, , F, D, F andfl, 

1 i 
s = s in — , c = cos — , 

2 2 

and the brackets designate the averaged values. Similarly, 

K,1 = (‘’I®') + 280s«c’)C'„,„ 

+ ^e’(60s’c=- 280s*c’)C;„„ 

*28s*)Ci,„ 

t28s’)C;„, 

+ ^ (1 + s' c2 (9 - 42 + 42 s’) C' 

.^e'sc3(l.7s'.fs’)c'„. 

+ I|e2s3c(10 -35 s2 + 28s^)C*212 

sc (1 - 9 s* + 21 s^ - 14 s®) 0^213^(1 - e2)-7/2 
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where 


and 


'“2010 ~ ^ '“20 ' 

^201 1 ~ ^20 ^ “ 2 . 

C'lii = +Cjj cos 2rr + S2j sin 2 tt , 

C'lij = + C'j cos { 2 tt - 4 Q) - Sjj sin { 2 tt - 4 Q), 
^2113 “ ■*■^21 2^ + 82 j sin 2 0 ., 

C'211 = - C'2 cos (2 rr - Q) 4 S'2 s in (2 77 - fi) , 
C2212 = + cos (277-3(7) 4 S'2 s in (2t7 - 3 ( 7 ) , 

C'213 = + Cj2 cos (7 - $22 sin ( 7 , 


[w;„] = 


+ |c2(i-ios2 + iss*)q,„, 


-|-s2 (6 - 20 s 2 4 ISs^*) c; 

O 


3102 


4 IsC (1 - 5 s2 4 Ss'*) C' 


3201 


^ 330 . 


4 — S^ C' 

8 3302 


+ ^scM- 3 s 2 )C',„, 


^3^c(2-3s2)q„j] e(l 


- 5/2 
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where 


^3101 =^3i '"o® + sinTf, 

^3102 =^31 ^o® ( 77 - 20 ) -S'j sin(7T- 2 0 ). 

^3201 “ ^32 S in (?7 - 0) , 

^3301 ■ ^33 (77 + 2 0 ) + S33 sin (tt + 20 ), 

^3302 = ^33 *^°® - S33 s in (77 - 40 ) . 

C3401 - C34 cos (77 + 0 ) - S34 s in (77 + 0 ) , 

C'402 = C34 cos (77 - 30 ) + S34 sin (77 - 3 0 ). 

[ W' j ] and [W'g] are trigonometric series in the arguments of the lunar 
theory and in 77 andO. The presence of 77 can be a source of resonances at the 
orbital inclinations previously mentioned. 

These resonances are of the same type as in the direct lunar effects and 
both effects cannot be separated from each other by means of harmomic analysis. 

We propose that in future planning the close neighborhood of the resonance 
zones be avoided. 

THE MAIN PROBLEM 

Making use of Lagrange equations: 

dSi _ 1 ^ 

n a* /l - sin i 

d&n _ 1 

n /l - e* s in i 
14 


and setting 


we obtain 


d Stt _ 1 - cos i dV »/l - 3V 

n /l - e^' sini^^ na^e 

dSM _ I - 3V 2 lv\ 

~ ” na^e *na^ \ / 


V = a^~ a'^ [W^q] 

nn 

= — a? a'^ k, x 

M 2 

(*-m ■) 

(‘ ‘I?'’) • 


dSi 

dt 



(3^) 
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■** ' *.(1 -e>)» “ 


n m 7 1 3 t 

• — a-* a k- X 


[-^( 

■-i-) 

cos i • C 

4 ( 

-ft-) 

ros i • C 

4 ( 

-S-) 

cos 2 i 
sin i 


<“ * (1 - e»)» “ 


n m 9 « 9 f 

• — a * 


(+3 / 

1 - ” 

rs \ 

42 

^3 / 

1 _ 3 2 


14 

^3 / 

1 - 

^2 \ 

; I?" 


(-1-2 cos i + 5cos^ i) C^qq 
( 1 - cos i) (3 + Scos i) CJjq 


I (*'17'’) (- * 

+ 2 • ^ a» a'» k, X 

at ** * 

*f (‘'if'*) 
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It is of interest to know the variation of the vector 


c = /l - R, 

and of the Laplace vector 

e = e P, 

where R is the unit vector normal to the osculating orbital plane, P the unit 
vector along the line of apsides directed toward the perigee, and Q = R x p 
completes the right-handed system. 

The differential equations of Milankovic (1939) for the variation of vectorial 
elements in our case become: 


d be 
dt 




3V _ 
— + e X 
Be 



d 5e 
dt 




BV 1 - e2 _ BV\ 

— + c X — 1 . 

Be / 


Taking into account that in the main problem 


S e = 0, 



the last system becomes: 


dSR _ 1 - BV 

n /l - 


( 4 ) 


d 8P 
dt 


1 


n a' 



( 5 ) 
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The mean values of the elements can be substituted into the right hand sides of 
these equations. Taking into account 


U| = f s in i s in n 
Uj = * 9 in i cos 0 

u^ = -t- cos i 
uj + u* + uj = 1 


we obtain from (3): 



Making use of 



— = i(I -RR) 
3c c ' 


( 6 ) 


where I is the idemfactor, we have 


and, because 


we have finally: 



+ 2V, 


3c 



Inserting this decomposition into (4) and (5), we obtain 


d8P 

dt 



1 

na^ /l - 



(7) 

(8) 


Substituting into (7) 


Rx -^ = 
3R 


/ Uj - Uj 3/3 Uj 

U 3 3/3uj - Uj 3/30^ 

YUj 3/3 Uj - Uj 3/3 Uj 


and taking (6) into account we deduce: 


d S Uj 

n 

m* 2 »3 1 _ 

• n X 

dt 

’ (1 - e»)» “ ’ 


fa/- 

55 ^ 2 \ 


W (■ 

"42 ^200 “2 ""3 



- (" ^21 ^2 ■*■ ^21 “3 


+ 1 (1 - H =') (‘'I - - Si, 
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d S u 

“d7 


2 n 

- = • — a k, X 

2n2 M 2 


{ 


(1 - e^V 

3 (. 55 

- -r (1 e 1 u, u, 

4 \ 42 / ^ ^ 


dS 


dt 


■-&•■) 

(+ Cjj Uj 

■-H-) 

[+ uj 

1 

n 2 
• — a'^ 

,2^2 M 

a'3 X 

1 - i- A 

14 1 

[2q, u, 

1 - ii A 

14 1 

(C 22 u, - 


or, in a form more convenient for integration, 
d 5 u. 




dt 


(1 -e^) 


2\2 M 


— k, X 


cos i2 s in i cos i 


Li fi .11, A c' 

\ 4 V 42 " j 
+ ^ ^ ^ ^ ^ ^21 ^ ■*■ ^21 


+ isin^ i (- Cj2 cos 2D. + S'j sin 2D) 
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d 8 u. 


n 2 »3 1 V- 

= + • — a-* a kj X 


dt (l-e2)2 M 


^ ^1 - s in i cos i (- Cji s in fl + Sjj cos H) 


d 8 u, 
dt 



55 

K " 42 


fi-i- 

^ 14 

-- 1 

2 ' 



(-1 


n 

(1 • 

-e2)2 

(.| 

fl- 5_, 
1 14 

1 

2 

fl-il, 

^ 14 


i cos i CjQQ s in n 


+ i s in^ i (+ C'j sin 2Q + Sjj cos 2Q) 


_ a2 a'3 k, X 
u 2 


sin^ i (- Cjj sin 2n + S'l cos 2fi) 
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Tables I-IV display the terms of the Fourier expansions of the tidal poten- 
tial, eq. (3), whose coefficients exceed 10“^. The effect of the time dependence 
of the obliquity, which must be taken into account for accurate long term pre- 
diction, manifests itself in the sixth decimal digit of the expansion coefficients. 

To obtain the total effect of the earth tides on an artificial satellite the 
differential equations of the Main Problem (3^) - (3*'^) must be supplemented 
by the indirect effect of the tides in the oblateness perturbations (Kozai, 1965) 
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TABI£ I. Development of C'200 
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TABLE II 
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TAfiU: m. Developoent of C'eao 
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d s n ^ ^ • /t • X 

=-i2tan i(di), 

dt 


d b 7T 

dt 


= (5 sin i 


-tan i)n(5i), 


dbM 

dt 


= 3n s in i(8 i) , 


where 81 is the tidal perturbation given by equation (3^) and (Brouwer, 1959) 


n = n 


r 3 JaR^cosi 
\ ^ a ^ (1 - 


3 R* cos i r , ^ 

+ ^ j^4 - 9 + 12 /l -e^ - (40 - 5 +36 /l - e^) cos^ ij 


32 a^ (1 -e^y 
15 (2+3e2) cos i 

32 a^(l-e2y 


(3-7 cos^ i) 


We have integrated the resulting set of differential equations with mean 
elements obtained for the BE-C satellite at the epoch, 1970 June 19.22547, and 
predictions for the tidal perturbations, neglecting the secular terms, are pre- 
sented in Figures I-IV. While an accurate determination of tidal parameters 
will require a longer span of data, comparison with available inclination resid- 
uals over a three month interval show to be in the neighborhood of .31. 

CONCLUSION 

The theory and the Fourier expansion of the tidal effects in the motion of 
an artificial satellite presented in this article permit a deeper Insight into the 
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nature of the tidal perturbations and allo^v a judgment about the relative signifi- 
cance of the amplitudes and periods of the terms in the series. 

The tidal effects are clearly noticeable in the observations and their com- 
parison with the theory will permit Improvement of the "global” Love numbers 
for the Earth. 

In our work we found inspiration in works by Kozai (1965) and Newton 
(1968) and also in work by D. Smith and R. Kolenklewlcz*, but we changed the 
method of expansion by resorting to the standard form of the lunar theory. In 

the present expansion all arguments are nearly linear in time. Thus the 
integration is easy and prediction for a longer interval of time is 
possible. 

We made use of the existing standard lunar theory because it is, presently, 
sufficient for our purposes. Any change in the coefficients, however, can be 
easily incorporated. 

The future addition to theory as the accumulation of observations and ex- 
perience proceed can go along several lines: more Love numbers can be in- 
cluded into the theory and the dependence of Love numbers on the longitude and 
latitude (Kaula, 1969) can be incorporated. 

ACKNOWLEDGMENT 

The authors wish to express their gratitude to their colleagues. Dr. D. 
Smith and R. Kolenklewlcz for interesting discussions and for sharing with us 
the results of their investigations and the data on the BE-C satellite. 

*private communication. 


27 



REFERENCES 


The American Ephemeris and Nautical Almanac, p. 539, Washington, 1972. 

Brouwer, D., Solution of the problem of artificial satellite theory without drag. 
Astron. Journ. 64, No. 9, pp. 378-397, 1959. 

Brown, E, W., Tables of the motion of the Moon, Yale University Press, 1919. 

Eckert, W. J., Walker, M. J. and Eckert, D., Transformations of the lunar co- 
ordinates and orbital parameters. Astr. Journ. 71, No. 5, pp. 314-332, 1966. 

Eckert, W. J., Van Flandern, T. C., Wilkins, G. A., A note on the evaluation of 
the latitude of the Moon. M. N. R. A. S. 146 . pp. 473-478. 

Fisher, D. and Felsentreger, T. L., Effects of the solar and lunar tides on the 
motion of an artificial earth satellite. Goddard Space Flight Center X-547- 
66-560, 1966. 

Kaula, W. M., Tidal friction with latitude dependent amplitude. Astr. Journ. 74, 
No. 9, pp. 1108-1114, 1969. 

Kozai, Y., Effects of the tidal deformation of the Earth on the motion of close 
Earth satellites. Publ. Astron. Soc. Japan, 17, p. 395, 1965. 

Milankovic, M., On the application of the vectorial elements in the computation 
of planetary perturbations (Serb.). Serb. Acad. Sc. M, 1939, Belgrade. 

Musen, P., Bailie, A. and Upton, E., Development of the lunar and solar pertur- 
bations in the motion of an artificial satellite, NASA TN D-494, 1961. 


Newton, R. R., A satellite determination of tidal parameters and earth decelera' 
tion. Geoph. J. R. Astr. Soc. 14, pp. 505-539, 1968. 

Siry, J. W., Astronomic and geodynamic parameters, NASA report 
X-550-70-481, 1970. 


29 









